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Introduction 


In this report, the progress we have made on the following projects during 
the period of this grant/ with, particular emphasis on the work carried out since 
submitting the semi-annual report in June i 

A. Calculation of the energy-density relation for pure neutron matter in the 

density range relevant for neutron stars using four different hard-core 
potentials. ■ 

B. Calculation of the properties of the superfluid state of the neutron component 
and the superconducting state of the proton component and the effects of 
polarization in neutron star matter. 



Neutron Matter - Normal State Calculations 

Since the preparation of the semi-annual status report, we have been able to 
obtain a better method of constrained minimization of energy. This new method gives 
energies considerably lower than those reported in the semi annual report. Therefore 
we have adopted this method and carried out the full range of calculations for Ohmura- 
Morita - Yamada-.(0MY-4) , standard hard core (SHC), Reid Hard Core (RHC) and Hamada- 
Johnston (HJ) potentials as outlined in the ■ proposal . 


The correlation factor used in our calculation has the following form; 

o) * } • ( c 

( o ' > r< 

The choice of this three -parameter correlation factor is based on considerable numerical 
experience involving the use of two simpler correlation factors; A two-parameter form 
(C2) which can be obtained from (C3) by choosing ,M( ^/^and a one-parameter form (Cl) 
obtained from (C 3 ) when - 0 . We have found that when we determine /' l { in 
by energy minimization, this is associated with the violation of several of the 
conditions discussed in the semi annual report, except in the low -density region. The 
additional flexibility of (C2) enables us to extend the region, where the conditions 
are satisfied, beyond the low densities. But (C2) turns out to be inadequate in the 
moderate to high density region. For example, it has been found in our earlier nuclear 
matter calculation using the SI1C potential that in order to extend the region of 
validity of the method up to and beyond the region of equilibrium density, we have to 
use the three -parameter form (C3). It should also be noted that the magnitude of the 
cluster correction £ relative to progressively diminishes as we moved from (Cl) to 

j & ■ 

(C2) to (C3) in our calculations. These remarks summarize some important general 

features of our numerical experience in neutron matter and symmetrical nuclear matter 



calculations . 
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The necessary conditions on the radial distribution function as well as other 
physically motivated conditions on the trial wave function discussed in the semi 
annual report are summarized below for convenience 


(lA) /} Kr- » n (Normalization condition £r structure-factor sum rule to 
lowest order) 

(ib) a 7\T ~ o 

(ll) 3l0)~ 0 ( Structure- factor sum rule) 

(ill) £(.(?.) ^0 (Structure-factor inequality) 

( IV ) 1<,£43 (Coulomb inequality) 

(VA) Xg - 0 (Pauli condition averaged over the fermi sea) 

(VB) r B - 0 (Pauli condition for the "average pair") 

We recall that some of these are mutually exclusive, if we attempt to impose all of 
them exactly on a state-independent correlation factor. However, we may try to choose 
a subset of these to be imposed exactly and hope to have the results thus obtained 
satisfy the remaining ones approximately. Our objective is to seek the best one 
among the various procedures that are possible within this general framework - best, 
in the sense that it gives the lowest energy. / 

In order to obtain a calculational method which is reliable in the high-density 
region, we chose 1 to try different procedures first at the density corresponding to 
kp ~ 3r5"Jm*using the OMY-4 potential. The difference in the methods is due to 
different choices of constraints on the energy minimization, as described below. 

Method Al. At each k^, the parameter ^ in the correlation factor A; r ) 

is determined by the condition ( IA ) , for every chosen set of values ( / } 

With - "'i thus determined, we find that for a given , energy £ — S fi ~t £■> t 
exhibits a minimum as a function of at some value The two-body approximation 

for energy, ( £_ also has a minimum with respect to y U f at approximately the same 

value /Uj . Therefore, we determine ~ / U ( by minimizing ( £p X £*) . The values of 
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} Y Z { Y) , S(o) and some other associated quantities defined in the semi- 
annual report, calculated for various values of ^ with z-3 are given in 
Table A-l. We then plot S(o) as a function of to determine where S(o)=0. 

But, it is seen from Table A~1 that S(o) never reaches zero. At somewhat lower 
densities , 3*0 S(o) crosses over from a negative value to a positive 


o 

value at some ^ . At we have to be satisfied with choosing K as the 

value where S(o) is closest to zero ; this value is = It* • At this , th 

other parameters fcY* and P, - Y* are redetermined from the condition (IA) and the 
minimization of £ respectively. In this manner, we arrive at the optimum set of 


parameters ( y* ^ $°) for the correlation factor at each density. For k F =3,SJY» ' 

and the 0MY-4 potential, the optimum set is ^ J 1 ^^7^) • the 

corresponding values for energy and other related quantities are entered under 
method (Al) in Table . We note that 0> which satisfies condition (IV), 

namely, the Coulomb inequality. By calculating &00 for a wide range of K values, 
we have also verified that condition (III), namely, the structure-factor inequality, 
is satisfied by the optimum correlation factor ^ j ° ^ also sa ''-'i- s ^ ^- eS 

the Pauli condition approximately since {X & (~ 6*0|l Finally, the correlation 


parameter It 1 = 0* U2 , is sufficiently small indicating good® convergence, especially 
since the ratios R c and R 5 (o} are also small. Method A2 . In this method, condi- 
tion (VA) is used Instead of (IA) to determine Y~Y otherwise the procedure in 
method Ai is followed. For the 0MX4 potential, the results obtained at Pp-’XSJXP 
using this method are given in Table A-2. We note that the energy obtained by this 
method, namely, £ = |2f^NeV£s higher than the value £~W)22 WeV obtained by method 
Al. Furthermore, this energy is associated with S(o)=-0.flSl compared to S(o)~ -0*007 
for method Al, Therefore method Al is better than method A2. 

Method A3 . The condition £=0 is used to determine Y^Y in this method, (In this 
case, 2. Aj\ r may be assumed to play the role of the "smallness parameter, ). Other— 
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wise the procedure in the two previous methods is followed* The results given in 
Table A-2 for lip- 3'5'fWt and .4 potential show that with E - |l5t?MtVand S(o)^ — 04-'2 / 

, this method is better than A2 , but still not as good as Al. 

We have carried out calculations using two other methods in which the con- 
dition determining % is changed to (IB) and (VB) respectively* These also re- 
sult in values of £ higher than that obtained in method Al and also associated 
with a greater degree of violation of the. condition S(o)=0* 

The method we adopted in calculating the results reported in this paper .is 
method A, which is an improvement over method Al, as we shall see below* In Al, 
at a given fcp, the value of which minimizes (£p + £2^ or £ for every chosen 
value of p is associated with S(o) < 0, except in the very low density region* 

However, at that stage of the procedure, we have ignored this violation of the 
structure-factor sum rule, anticipating that \ J 2 will be adjusted to make S(o)=0, 

In studying S(o) as a- function of p for given p, , we have observed the behavior 
schematically indicated in figure A-l, The five sketches are for five different 
values ( of namely j3 , such that * ^ ur nujllerica - j - results show 

that the region near the first .(deeper) minimum in S(o) is associated with ex- 
tremely large values of "g and severe violations of all the conditions, except the 
one explicitly imposed to determine ~i * We are seeking the lowest value of 
consistent with S(o)=o» But crisis not an acceptable value of p even though S(o) 
is zero there, because of the reasons mentioned above* As the curves (a), through 
(e) in figure A-l indicate, the value p of ^ where £ is minimum is roughly 

1 ■ 

the same location where S(o) also has a minimum, which is negative unless is 
sufficiently large* In method Al, we take a high enough value of to make 
S(o)*=0, so that the situation in curve (d) is obtained* In general, the 'energy £ 
increases with . In order to see whether we .can get a lower value of £ than 
obtained in method Al, we attempt to determine p by not merely minimizing 
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but by seeking its minimum consistent with S(o)=o. For the value , we see that 

this is not possible because ('!) we have ..already seen that c£ ( is not acceptable and 

V 

( 2 ) the only other zero of S(o) occurs at a very high pj where, once again, the 
value of is found to be too large and some of the conditions are violated. ^It 
may be appropriate here. to point out that we have found the following general result 
in our numerical calculations using various methods: the coulomb inequality and the 

structure-factor inequality along with the conditions that]j,'Xg and A if be small would 
restrict the parameter space to a region where the difference between the values of 
(-{ and is not very large. Thus we find, referring to curves (a) through (e) , 
that the zero of S(o) occurring at a very high V 


and the zero of S(o) at a Very small f~s - ot f are both unacceptable.^' 


Thin 


we are forced to go to higher values of . For p^ « p, (curve (b), S(o) is zero 

at p ( — o( 0 ♦ Note that since £(p ( ) is the minimum of £()->), Now, if we 

go to curve (c) with S(o) is zero at two values of , namely and . The 

energies are related as £ (o( 2 ) >£6^) > £ ( p t ) . In curve (d), for j-i, -■ l%> S(o) - o 

when is and ot^ with the relation between the energies > £ (o£ 2 )>£ For 

values of y greater than , S(o) has only the zero 'at y > Since the second 

r 1 

minimum in S(o) is already on the positive side and it becomes more positive with 
increasing . This situation is illustrated in curve (e). 

Thus, in method A, for each p^ , we are determining p' by • minimizing £ subject 
to the condition, S(o) =0. The remaining task is to find the minimum of £ in the 
range of values spanned by £(a(c)f £ (‘^3) and £ (dq) which are all consistent with the 
structure- factor sum rule. At some densities we find that the minimum of occurs 

when Pp, ot c • At others, the minimum of £ (p^) . At still other densities the 

value of ^ at which is lowest lies between o( o and oC^. Denoting this minimizing 

value by y® , we have ( p/ c ; for the optimum set of parameters, where is 

the value of at whi.ch has the lowest value consistent with S(o)=o, as 
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explained above. The results for £ ( ^ j ^ ) an ^ t ^ ie associated quantities at 

k„ = 3-5 fm" 1 obtained using the OMY-4 potential are given in table A-2. We note 

r 

the significant improvement in the results over those from method Al, especially in 

the value of energy, namely, £. = 722*5 Mev compared to £■ = 1022 Mev from method Al. 

The results from the calculation for the OMY-4 potential using method A for the 

entire density range 0.25 f m kp 5*00 fm ^ are given in table A-3* The values 

of £ determined by this method are plotted against the density parameter in figure 

A-2. The most notable feature of the£(k ) curve is the occurrence of a local mini- 

r 

mum near k . = 1-75 Tm • ■ - > 

r 

and the discontinuity in the slope of £(l^) near = 1.30 fm . The energies at all 
densities, except in the very low density region, determined by our constrained mini- 
mization procedure have come out much higher than the results obtained by other authors. 
The convergence of the cluster series for energy and other associated quantities appears 
to be excellent judging from the data in table A-J. One of the most important Aspects 
of our method that emerge from these results is the crucial role played by the structure 
factor sum rule in achieving this convergence. It is also responsible forq-aising the 
energies to such high values. For example, if the condition s(o) = o is ignored, at 
kp = 1.5 fm ' , we will get £ = 13 Mev compared to £ = 23 Mev we have in table A-3, for 
the OMY-4 potential. This effect, of course, becomes more pronounced at higher 
densities. Thus the answer to the question whether the energies given by method A are 
reliable depends crucially on the justification for the central importance we have 
assigned to the structure- factor sum rule. (We may- also point out here that once S(o) 

=o is imposed in the manner we have done in method A, S(k)>_ o is always automatically 

satisfied. Choosing one of the other zeroes of S(o) leads to a violation of S(k)>>Q 
for some range of k values and, associated with this,, to large values of ^ , except in 

the very low density region, where a straightforward minimization of Mj) with 

respect .to; yU.^ and yields a minimum that is consistent with all the conditions.) 

Since our method is variational, it may be contended that 



7 


our high values of energy may constitute only a poor upper bound to the true energy * 

In reply to this criticism, we can only say that our numerical experience as summarized 
in this report suggests that this is the lowest energy we can obtain for- our choice 
of constraints with the three-parameter correlation factor (CT5). We also note that 
the constraints chosen in method A ;are not arbitrary, but have been justified by the 
physical arguments referred to in the original proposal and by the numerical evidence 


for their central role in obtaining apparently rapid convergence of the energy series 
and the associated cluster expansions. 

The results from the various methods we have attempted indicate that any one of 
the conditions (lA), (IB), (VA) and (VB) tend to make the correlation parameter 
small. But we also find that smallness of g by itself does not assure the convergence 
of the cluster expansions. Only if the smallness of is accompanied by the satis- 
faction of S(lt)> o, it seems possible to obtain over-all convergnece. In spite of 
this statement, g still plays a useful role as an ordering parameter for cluster 
expansions in the CBF formalism. We have repeated our calculation at =3*5 fm”^ 
for OMY-4 potential using method Al, omitting the term 8^ from the expression for 8^ - 
The energy obtained in this way is higher than the value £ = 1022 Mev given, in table 
A-2. We recall that £ is a (reducible) four-body correlation term whereas all 


other terms in contain the effects of (irreducible) three-body correlations. If 
we are to use the "number of bodies" to classify various orders in the cluster series, 
we must exclude from <£ ^ . But, according to the - classification ' scheme, 
is of the same order as the other terms constituting and hence must be included in 
e Hence the numerical results mentioned above seem to justify the use of as an 
ordering parameter for cluster expansions. 

The energy per particle of neutron matter calculated using the SHC, RHC and HJ 
potentials are given in figures A- 3 , A -4, and A-5 respectively. These also show local 
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minima and discontinuity in slope as we have seen for the OMY-4 potential. We 
cannot yet. attach any physical significance to these nonmonotonic regions of the 
curve for the following reason. The parameter values in the correlation factor 
vary significantly- .with density in these regions as seen from table A~3 for the 
case of the OMY-4 potential. It is possible that when perturbation corrections 
are calculated using these optimum correlation factors, the final result, namely, 

(£ f +£ 7 ) +&Z when plotted against kp, may not have such local minima or 

discontinuities in its slope. When is calculated on the basis of a 

convergent perturbation expansion, if it is found that the final energy- density 
curve still displays these features, then only we can consider these as possible 
indications of phase transitions. 

In the meantime, we also need to look critically at our choice of the structure 
factor-sum rule, namely S(o) = o, as a vital constraint in the energy minimization. 

In particular, we need to assure ourselves that it is realistic to impose that 
condition on a short range correlation factor and that the raising of energies caused 
by it is not an artificial effect. We want to emphasize, however, that we have given 
some physical arguments that justify the use of -this condition and also the fact that 
it definitely helps to obtain better convergence of the cluster expansions in our 


formalism. 
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B „ Effects of Polarization on Neutron Star Matter 

The final numerical work on this preliminary study of this subject has been 
completed. We plan to present a paper on this subject at the APS meeting to be held 

in Washington, D. C. next April. 

We have completed the following additional calculations since our last semi- 
annual report: 

1. The three-particle contributions to the normal state energy (per particle) 
in the density region interested to us has been carefully checked out. 

They are as following: 


At density (expressed 
in term of (fm - 1 )) 

0.25 

0.50 

1.00 


the three particle 
‘ contributions, ^(MeV) 

0.0000 

- 0,0012 

-O.O858 


These numerical results for normal state indicate that the three-particle con- 
tributions in the density region interested to us are indeed sufficient small as 
compared to their corresponding single and .two particle contributions to the energy ex- 
pectation value. That leads us to believe that the three-particle contributions 
to the energy expectation value is ngligible as compared with the smgle-and two- 
particle contributions in the superfluid state because of the same cluster ex- 
pansion employed in both cases. 

b. The numerical accuracy of the condensation energy calculations were tested by 
.-changing our computer programs for those calculation from the 12-point Gaussian 
integral to the 20-point Gaussian integral for the enhanced factor £ = 1.50 for 
the calculations mentioned in our last annual report (see Table B-2) . 
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The results are as following: 


at density, k^, (fm”^) 


Sc CMeV) 


12- PT.. G. I. 

20-pt. G. I. 

Difference 

0.60 

2 .O 695 

2.0642 

0.0053 

O .96 

1 . 1047 

1.1095 

0.0048 

Thus the numerical evaluations 

for condensation 

energies are reliable 

up to the 


third place after decimal point (the difference is less than one -half of one percent) 
3. In addition to the numerical results given in the last two reports, we have 
made the following two sets of calculations to complete our preliminary study 
of this subject: 

(a) The condensation energy is calculated in such a way that it has the effect 
of incorporating polarization prior to the short range correlations and 
superfluidity. Numerically, the condensation energy, £ c , is calculated 

by using the optimal energy gap, 4 k , f° r each given enhanced ' factor f6 

.. in both £ c and £ 3 but only the in Pkfc is multiplied by A . The 
{ c 2/ ’ 
results are given in Table (B-l) for enhanced factor ^3 = 1-20, 1.30 and 

1.40. 

(b) The condensation energy is calculated in such a way that it has the 

effect of incorporating the short range correlations and superfluidity 
prior to polarization. Numerically, the condensation energy, , is 
calculated by = £ c + (i.e. the whole <^in P k ^ is enhanced by the 

factor p). The results are given in Table (B-2) for enhanced factor 

1.15, 1.20, 1.30 and 1.40. 

The results of Table (B-l) of this report and Table (B-2) of the last annual report 
all seem to indicate that the polarization effect indeed enhances the condensation . 
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energy and the gap function and there is a tendency of the neutron star matter to 
undergo a first-order phase transition at a relatively lower density region 
(10 12 -10 1 ^ gm-cm" 3 ). This preliminary study reaffirms the necessity and interest 
of a first principle theory study of the problem. 



Table (B-l) <J 

In this set of calculations, the condensation energy, £ , is calculated by using the (^fej . . for each corresponding Bin both £ c and £ c , but only the ft 0 \vu 
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